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has been found. We investigate the role that the B* and Bq{Q^) resonances might play in 
these analyses. 
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The decay mode B ^ pir is currently studied as a channel allowing, in principle, to measure 
without ambiguities the angle a of the unitarity triangle. It is also investigated by the CLEO 
Collaboration where a branching ratio larger than expected for the decay mode B^ — » p'^Tr* has 
been found. We investigate the role that the B* and _Bo(0+) resonances might play in these analyses. 



I. INTRODUCTION 



The measurement of the angle a in the unitarity triangle will be one of the paramount tasks of the future b-factories, 
such as the dedicated e+e" machines for the BaBar experiment at SLAG and the BELLE experiment at KEK 
or hadron machines such as the LHG at CERN, with its program for B-physics Q Differently from the investigation 
of the /? angle, for which the B — > J/'^Ks channel has been pinned up Q and ambiguities can be resolved the task 
of determining the angle a is complicated by the problem of separating two different weak hadronic matrix elements, 
each carrying its own weak phase. The evaluation of these contributions, referred to in the literature as the tree 
(T) and the penguin (P) contributions, suffers from the common theoretical uncertainties related to the estimate of 
composite four-quark operators between hadronic states. For these estimates, only approximate schemes, such as the 
factorization approximation, exist at the moment, and for this reason several ingenuous schemes have been devised, 
trying to disentangle T and P contributions. In general one tries to exploit the fact that in the P amplitudes only the 
isospin-1/2 Qpart of the non-leptonic Hamiltonian is active |^; by a complex measurement involving several different 
isospin amplitudes, one should be able to separate the two amplitudes and to get rid of the ambiguities arising from 
the ill-known penguin matrix elements. 

One of the favorite proposals involves the study of the reaction B — > pi:, i.e. six channels arising from the neutral 
B decay: 

B° ^ p+n- , (1) 
B" ^ p-7r+ , (2) 
B° - , (3) 

together with the three charge-conjugate channels, and the charged decay modes: 

B- ^ p-TrO , (4) 
B- ^ , (5) 

with two other charge-conjugate channels. Different strategies have been proposed to extract the angle a, either 
involving all the decay modes oi a. B into a pir pair as well as three time-asymmetric quantities measurable in the 
three channels for neutral B decays |^-|^, or attempting to measure only the neutral B decay modes by looking at 
the time-dependent asymmetries in different regions of the Dalitz plot^. 



^Opportunities for _B-physics at the LHC have been recently discussed at the workshop on Standard Model 

Physics (and More) at the LHC, 14-15 October 1999; copies of transparencies can be found at the site 

http://home.cern.ch~mlm/lhc99/octl4ag.html. 

— n 

If one neglects electroweak penguins. 

^In this way the measurement of a decay mode with two neutral pious in the final state, eq. (0), can be avoided. 



1 



Preliminary to these analyses is the assumption that, using cuts in the three invariant masses for the pion pairs, 
one can extract the p contribution without significant background contaminations. The p has spin 1, the tt spin as 
well as the initial B, and therefore the p has angular distribution cos^ 9 {0 is the angle of one of the p decay products 
with the other tt in the p rest frame). This means that the Dalitz plot is mainly populated at the border, especially 
the corners, by this decay. Only very few events should be lost by excluding the interior of the Dalitz plot, which is 
considered a good way to exclude or at least reduce backgrounds. Analyses following these hypotheses were performed 
by the BaBar working groups MonteCarlo simulations, including the background from the /o resonance, show 
that, with cuts at m^^ = rUp ± 300 MeV, no significant contributions from other sources are obtained. Also the role 
of excited resonances such as the p' and the non-resonant background has been discussed Q . 

A signal of possible difficulties for this strategy arises from new results from the CLEO Collaboration recently 
reported at the DPF99 and APS99 Conferences 

B{B^ p°7r±) = (1.5±0.5±0.4) X 10"^ , (6) 

B{B p^TT*) = (3.5l};J±G.5) X IQ-s , (7) 

with a ratio 

_ B{B ^ p^TT^) _ 2 3+1 3 (f^) 

^ - B{B± - pOn±) ^ 2.3 ±1.3. (8) 

As discussed in [Q, this ratio looks rather small; as a matter of fact, when computed in simple approximation schemes, 
including factorization with no penguins, one gets, from the DDGN model of Rcf. R ~ 13, admittedly with a 
large uncertainty; another popular approach, i.e. the WBS model ||l2|, gives i? ~ 6 (in both cases we use ai — 1.02, 
02 = 0.14). The aim of the present study is to show that a new contribution, not discussed before, is indeed relevant 
to the decay (||) and to a lesser extent to the decay (||) . It arises from the virtual resonant production depicted in Fig. 
1, where the intermediate particle is the B* meson resonance or other excited states. The B* resonance, because of 
phase-space limitations, cannot be produced on the mass shell. Nonetheless the B* contribution might be important, 
owing to its almost degeneracy in mass with the B meson; therefore its tail may produce sizeable effects in some of 
the decays of B into light particles, also because it is known theoretically that the strong coupling constant between 
B, B* and a pion is large IT^ . Concerning other states, we expect their role to decrease with their mass, since there 
is no enhancement from the virtual particle propagator; we shall only consider the 0+ state Bq with = 0+ because 
its coupling to a pion and the meson B is known theoretically to be uniformly (in momenta) large [l^ . The plan 
of the paper is as follows. In Section || wc list the hadronic quantities that are needed for the computation of the 



widths; in Section III we present the results and finally, in Section IV, we give our conclusions. 



II. MATRIX ELEMENTS 

The effective weak non-leptonic Hamiltonian for the |Ai3| = 1 transition is 



(9) 



The operators relevant to the present analysis are the so-called current-current operators: 



{daUp)v-A{ui3ba)v-A 

{daUa)v-Aiupbi3)v-A , (10) 



and the QCD penguin operators: 



*We omit, as usual in these analyses, the electroweak operators Qk (k — 7, 8, 9, 10); they are in general small, but for Qg, 
whose role might be sizeable; its inclusion in the present calculations would be straightforward. 
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Qs = {daba)v-A ^ {q'pq'p)v-A 

q' —u,d.s,c,b 

Qa = idabi3)v-A ^ {q'pq'Jv~A 
q' —u.d.s.c.h 

Qs = {dc,ha)v-A ^ {q'i^q'i^)yj^A 

q' —u,d.s,c,b 
q' —u.d.s.c.b 



(11) 



1.100, C3 = 0.012, C4 = -0.029, 
: 225 McV, n = fhb(mb) ^ 4.40 GcV 



We use the following values of the Wilson coefficients: Ci = —0.226, C2 = 
C5 = 0.009, Ce = -0.033; they are obtained in the HV scheme @, with A^|^^ 

and mt = 170 GeV. For the CKM mixing matrix |l5| we use the Wolfenstein parameterization with p = 0.05, 
•q = 0.36 and A = 0.806 in the approximation accurate to order A'^ in the real part and in the imaginary part, i.e. 
Vud = 1- AV2, Kb = AA3 [p-^77 (1 - AV2)], Vtd = AX'^ {1 - p - ii^) and Vtb = 1. 

The diagram of Fig. 1 describes two processes. For the B* intermediate state there is an emission of a pion by 
strong interactions, followed by the weak decay of the virtual B* into two pions; for the p intermediate state there is 
a weak decay oi B ^ pn followed by the strong decay of the p resonance. We compute these diagrams as Feynman 
graphs of an effective theory within the factorization approximation, using information from the effective Lagrangian 
for heavy and light mesons and form factors for the couplings to the weak currents ^. 

To start with we consider the strong coupling constants. They are defined as 



{B^p') 7r-iq)\B*-ip,e)) 
{B-{p') I 
(^"(q')^-(<z) I 



ffpe ■ {q' - q)- 



(12) 



In the heavy quark mass limit one has 



9 



mBoms s 



(13) 
(14) 



For g and h we have limited experimental information and we have to use some theoretical inputs. For g and h 
reasonable ranges of values are g = 0.3-0.6, —h ~ 0.4-0.7 [|l^. These numerical estimates encompass results obtained 
by different methods: QCD sum rules ||l3|, potential models jlj], effective Lagrangian ||l^, NJL-inspired models [ pO[ . 
Moreover gp = 5.8 and ~ 130 MeV. This value of gp is commonly used in the chiral effective theories including the 
light vector meson resonances and corresponds to V p ~ 150 MeV; see, for instance, where a review of different 
methods for the determination of g is also given. 

For the matrix elements of quark bilinears between hadronic states, we use the following matrix elements: 



2ma 



(tt- |d75u| 0) 
{n'iq) \u^rM B* - {p)) = ze^^ {q - p) , 



{p+{q, e) \u^,b\ &{p))=ze*^{p~q), 



2mB'Al 
rrib + niq 
2mpAo 



rub + ruq 



{TT+iq) \u^pb\ B°{p))=F, 



[p + qT- 



— m. 

{p - q)' 



Hp-qr 



Fn 



{p - q) 



rip-qY 



^In the second reference of Q a similar approach has been used to describe the decay mode B'^ — > D^D~-k'^; the main 
difference is that for B — > Stt we cannot use soft pion theorems and chiral perturbation theory, because the pions are in general 
hard; therefore we have to use information embodied in the semileptonic beauty meson form factors. This is also the main 
difference with respect to M. 
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(7r+(g) |S7^(1-75)6| = * 

(0 \urd\ p-{q, e)>=/pe^ 



'Bo 



iv-q) 



~ m 
Fo 



^Hp-qY 



where /p = 0.15 GeV^ 



and 



^0 



Al{Q) = 0.16, 
Fi=Fi(0) = Fo(0)-0.37 



Ao = ^o(O) = 0.29 , 
Fo- = Fi(0)-Fo(0) 



-0.19 



(15) 



(16) 
(17) 



The first three numerical inputs have been obtained by the rclativistic potential model; and Fi can be found in 
pl| , while Aq has been obtained here for the first time, using the same methods. The last figure in ( [l7| ) concerns 
Fq , for which such an information is not available; for it we used the methods of ||l^ and the strong coupling BBqit 
computed in |23]. 



III. AMPLITUDES AND NUMERICAL RESULTS 



For all the channels we consider three different contributions Ap, Ab-', Abq, due respectively to the p resonance, 
the B* pole and the _Bo positive parity 0+ resonance, whose mass we take ^ to be 5697 MeV. 
For each of the amplitudes 



A--+ = A{B- TT~TT~TT+) 
= A{B- ^ TT-Tr°TT°) 

= A(B" ^ ^+7r-7r") 



(18) 
(19) 
(20) 



we write the general formula]] A*-'* 



P 



A 



ijk 



A^^ . We get, for the process (^Sl 



A — + 



t — u 



t ~ m?p-{- iTpnip t' — m'^ + iTpTrip 

U{t,u) n{t',u) 



A 



BO 



t — in^gt + iT B'''niB' t' — m^, + iT b'Ttib* 



* ~ "^Bo + ^^B„mBo t' - m^^ + iVBoinBa 



(21) 



where, if p,r- is the momentum of one of the two negatively charged pions t = (p^- +p^+)^, t' is obtained by 
exchanging the two identical pions and u is the invariant mass of the two identical negatively charged pions. Clearly 
one has u + t + t' ~ m'^g + 3m^. The expressions entering in the previous formulas are 



71 "'71 



fpFi[ci + ^] +?7ipAo 



Gp 

71 



VtbV, 



*9p_ 
V2 



td 



(c4 



C3 



ifpFi 



pAoA) + 2 (c6 + I) mpA,f ^""'.^ 

V 3 / (mb + mqjimq 



K = 



if" 



71 

mBo " V2 



aV2 g mlAZ {KbKTd (^2 + y ) - ^t^K 

4) {VubV:^ (c2 + y) - Vt,V^ 



3 ( ^ ^ 3 V (mf, + mq)2mq 



rriB " ^ / 



m: 



C3 _ 2 + ^ii") .. 

3 V 3 / (jTif, + mq)2mq 



} , (22) 



with mt, = 4.6 GeV, m„ 



6 MeV, Fs* = 0.2 keV, Fb„ = 0.36 GeV [0. Moreover, for the process (|T|) 



® We identify the 0^ state mass with the average mass of the B** states given in [g4|. 

'^We add coherently the three contributions; the relative sign the B resonances on one side and the p contribution on the other 
is irrelevant, as the former are dominantly real and the latter is dominantly imaginary. The relative sign between B* and Bq 
is fixed by the effective Lagrangian for heavy mesons. 
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A 



A 



-00 _ 

p 


V 


-00 


1 


B' — 




-00 




Bo - 





1 r s + s'-Aml KIl{s\s) + /vin(s',u) KIl{s,s') + A'in(s,u) 
— < K 1 = 1 = 4- 



2m|. 



s' — m^, + iTB'niB* 



s — m^, + iVs'TTiB 



- ml) 



^Bo +i^BomBo 



m 



B„ +i^BomBo 



4n + ^^B„mB„ 



(23) 



In this case we define s ~ {p^- 



if Pt^o is the momentum of one of the two identical neutral pions, s' is obtained 



by exchanging the two neutral pions and u is their invariant mass (again we have a relation among the different 
Mandelstam variables: s + s' + u = m\ + 3to^). Then 77", Ki and K'^'^ are given by 

Gp 



1] 

and 



^/2 
71 



* 9p 
9p 



Cl 



mpAof^ Cl 



C2 



3 / 



{-fpFi + mpAafrr) - 2 ce + — mpAof^ 



C5 



C4 



C3 



2 C6 



C2 



V2^ 



C4 



C5 

' 3 
C3 



(to;, + mq)2mq 

™2 



(mb + mq)2mq 
-2 fee 



C5 
3 



(mt + mq)2mq 



} , (24) 



n(.T,2;) = - - 

while k'^ was given above in (p2|). 

Finally, for the neutral S decay (EQ), we have 



V , - - x) 



(25) 



t — u 



t - 7712 + iFp mrip 



/I+-0 



-Do 



A'n(s,f) + Kili{s,u) 

s — m?g, + iTB'niB* 

( k° + A''^^ 



t — m-^, + iYB'TRB 



■5 - '"so + ^^Bo-niBo t - mi + iTBarriBo 



{ml -ml) , 



(26) 



where s = (p^- +PTro)'^, t = {p^- + p-^^)"^, u = (pjr+ +?57ro)^i and s + t + u^ m\ + 377i^. The constants appearing in 
these equations are : 



= -^ifpFi + rUpAoU)^ [VutV:, (ci + |) + Vt,V;, (c4 + |) 



+ 



G 



V2 



C5 



3 / (jTib + mq)2mq 



77^ 



G \ c 
gpfpF,^ {k^K*. (c2 + y) - ^tfc^;. (c4 + y) } 



C4 + y - 2 ( C6 



(mb + mq)2mq 



(27) 



For the charged B decays we obtain the results in Table I and 11. In order to show the dependence of the results 
on the numerical values of the different input parameters, we consider in Table I results obtained with g = 0.40 
and h = —0.54, which lie in the middle of the allowed ranges, while in Table II we present the results obtained 
with g = 0.60 and h ~ —0.70, which represent in a sense an extreme case (we do not consider the dependence on 
other numerical inputs, e.g. form factors, which can introduce further theoretical uncertainty). In both cases the 
branching ratios are obtained with tb = 1-6 psec and, by integration over a limited section of the Dalitz plot, defined 
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as nip — 5 < (Vi, \/t) < nip + 5 for B~ — ^ 7r~7r~7r+ and nip — 5 < {^/s, < nip + S for B~ — > 7r~7r°7r°. For 5 
we take 300 MeV. This amounts to require that two of the three pions (those corresponding to the charge of the p) 
reconstruct the p mass within an interval of 26. Numerical uncertainty due to the integration procedure is ±5%. 

We can notice that the inclusion of the new diagrams {B resonances in Fig. 1) produces practically no effect for the 
B^ — > 7r^7r°7r° decay mode, while for B~ tt~^tt~tt~ the effect is significant. For the choice of parameters in Table 
I the overall effect is an increase of 50% of the branching ratio as compared to the result obtained by the p resonance 
alone. In the case of Table II we obtain an even larger result, i.e. a total branching ratio B{B~ tt^tt~tt~) 
of 0.82 X 10~^, in reasonable agreement with the experimental result (||) (the contribution of the p alone would 
produce a result smaller by a factor of 2). It should be observed that the events arising from the B resonances 
diagrams represent an irreducible background, as one can see from the sample Dalitz plot depicted in Fig. 2 for the 
B~ 7r+7r~7r~ (on the axis the two _ squared invariant masses). The contributions from the B resonances 

populate the whole Dalitz plot and, therefore, cutting around t ^ t' ^ nip significantly reduces them. Nevertheless 
their effect can survive the experimental cuts, since there will be enough data at the corners, where the contribution 
from the p dominates. Integrating on the whole Dalitz plot, with no cuts and including all contributions, gives: 

Br(B" ^ 7r"7r"7r°) 1.5 X 10"^ Br(B" ^ tt+tt^tt") = 1.4 x 10"^ (28) 

where the values of the coupling constants are as in Table I. 

We now turn to the neutral B decay modes. We define effective width integrating the Dalitz plot only in a region 
around the p resonance: 

T^ffiB" ^ p-n+) = r(B° ^ n+n-n") (29) 
T,ff{B°^p+n-) = r{B"^n+n-7T") , (30) 

T^niB" - P°7r°) = r(B" ^ n+n-n') (31) 

The Mandelstam variables have been defined above and again we use i5 = 300 McV Similar definitions hold 

for the B^ decay modes. The results in Table III show basically no effect for the B^ p^n^ decay channels and 
a moderate effect for the p^Tr" decay channel. The effect in this channel is of the order of 20% (resp. 50%) for _B" 
(resp. i?") decay, for the choice g = 0.60, h = —0.70; for smaller values of the strong coupling constants the effect is 
reduced. Integration on the whole Dalitz plot, including all contributions, gives 

Br(i?° ^ n+TT-n") = 2.6 x 10~^ (32) 

confirming again that most of the branching ratio is due to the p-exchange (the first three lines of the p column in 
Table III sum up to 2.3 x lO"'^). 

To allow the measurement of a, the experimental programmes will consider the asymmetries arising from the 
time-dependent amplitude: 

Ait) = e-^/'' (cos^A+-° ±2sin^A+-A , (33) 



where one chooses the ± sign according to the flavor of the B, and Am is the mass difference between the two mass 
eigenstates in the neutral B system. Here A"' ° is the charge-conjugate amplitude. We have performed asymmetric 
integrations over the Dalitz plot for three variables: i?2 and i?3, which multiply, in the time-dependent asymmetry, 
respectively 1, cos Amt and sin Amt. We have found no significant effect due to the B* or the Bq resonance for Ri and 
i?3. On the other hand these effects are present in i?2, but i?2 is likely to be too small to be accurately measurable. 



IV. CONCLUSIONS 



In conclusion our analysis shows that the effect of including B resonance polar diagrams is significant for the 
B^ — > 7r~7r~7r"'" and negligible for the other charged B decay mode. This result is of some help in explaining the 
recent results from the CLEO Collaboration, since we obtain 

R = 3.5 ±0.8 , (34) 
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to be compared with the experimental result in eq. (|[) . The p resonance alone would produce a result up to a factor 
of 2 higher. Therefore we conclude that the polar diagrams examined in this paper are certainly relevant in the study 
of the charged B decay into three pions. 

In the case of neutral B decays we have found that, as far as the branching ratios are concerned, the only decay mode 
where the contribution from the fake p's (production of a pion and the B* or the Bq resonance) may be significant is 
the neutral p^ir^ decay channel. As for the time-dependent asymmetry no significant effect is found. Therefore the 
B — !■ TTTTTT decay channel allows an unambiguous measurement of a, with two provisos: 1) only the neutral B decay 
modes are considered; 2) the p^Tr" final state can be disregarded from the analysis. 
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Fig. 1 - The polar diagram. For tlic B resonances {B* = 
1~ , 0^) the strong coupling is on the left and the weak 
coupling on the right; the situation is reversed for the p 
production. 



TABLE III. Effective branching ratios for the neutral B 
decay channels into pvr {g — 0.60, h — —0.70). Cuts as indi- 
cated in the text. 



Channels 


P 




p + B* 




p + B* + Bo 




0.50 X 10 


— 5 


0.52 X 10 


— 5 


0.49 X 10"^ 


B° ^ P+-K- 


1.7 X 10" 


5 


1.7 X 10" 


5 


1.7 X lO"'^ 


B" p\" 


0.10 X 10 


-5 


0.15 X 10 


— 5 


0.12 X 10"^ 


-> p+TT- 


0.49 X 10 


-5 


0.51 X 10 


— i) 


0.48 X 10"" 


B° ^ p-7r+ 


1.7 X 10~ 


5 


1.7 X 10" 


5 


1.7 X 10"^ 


B" ^ p%° 


0.11 X 10 


— 5 


0.17 X 10 


— 5 


0.15 X 10"^ 




Fig. 2 - Sample Dalitz plot for the decay B^ 
TT+TT^TT^. In order to show the mass distribution of the 
B resonance diagrams, only their contribution is taken 
into account for this plot. 



TABLE L Effective branching ratios for the charged B de- 
cay channels into three pions for the choice of the strong cou- 
pUng constants g = 0.40 and h = —0.54. Cuts as indicated in 
the text. 

Channels p p + B* p + B* + Bo 

B- -> 7r-7rV° 1.0 X 10"^ 1.0 x 10"^ 1.0 x lO"'"^ 
B- ^ 7r+7r"7r" 0.41 x 10"^ 0.58 x 10"^ 0.63 x 10"^ 



TABLE II. Effective branching ratios for the charged B 
decay channels into three pions for the choice of the strong 
coupling constants g = 0.60 and h — —0.70. Cuts as indicated 
in the text. 



Channels 


P 


p + B* 


p + B* + Bo 


B- 7r-7rV° 


1.1 X 10"^ 


1.0 X 10"^ 


1.1 X lO-'^ 


B" tt^tt^tt" 


0.41 X 10"^ 


0.74 X lO"'^ 


0.82 X 10~^ 
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